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ABSTRACT 

We investigate the 21 cm signals from the UV ionizing sources in the reion- 
O ■ ization epoch. The formation and evolution of 21 cm emission and absorption 

regions depend essentially on the kinetics of photons in the physical and frequency 
spaces. To solve the radiative transfer equation, we use the WENO algorithm, 
which is effective to capture the sharp ionization profile and the cut-off at the 
front of light (r = ct) and to handle the small fraction of neutral hydrogen and 
helium in the ionized sphere. We show that a spherical shell of 21 cm emission 
and absorption will develop around a point source once the speed of the ion- 
ization front (I-front) is significantly lower than the speed of light. The 21 cm 
shell extends from the I-front to the front of light; its inner part is the emission 
^ ■ region and its outer part is the absorption region. The 21 cm emission region 

depends strongly on the intensity, frequency-spectrum and life-time of the UV 
ionizing source. At redshift 1 + z = 20, for a UV ionizing source with an intensity 
E ~ lO^ergs" 1 and a power law spectrum v~ a with a = 2, the emission region 
has a comoving size of 1 - 3 Mpc at time ~ 2 Myr. Nevertheless, the emission re- 
gions are very small, and would be erased by thermal broadening if the intensity 
is less than E ~ 10 43 ergs _1 , the frequency spectrum is thermal at temperature 
T ~ 10 5 K, or the frequency spectrum is a power law with a > 3. On the other 
hand, the 21 cm absorption regions are developed in all these cases. For a source 
of short life-time, no 21 cm emission region can be formed if the source dies out 
before the I-front speed is significantly lower than the speed of light. Yet, a 21 
cm absorption region can form and develop even after the emission of the source 
ceases. 

Subject headings: cosmology: theory - intergalactic medium - radiation transfer 
- methods: numerical 
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1. Introduction 

The detection of redshifted 21 cm signals from the early universe is attracting many 
attentions in the study of cosmology, because it provides a window to probe the baryonic 
gas and the first generation of light sources in the cosmic dark ages (e.g. Furlanetto et al. 
2006). Observational projects of highly sensitive meter radio telescopes, such as the Low 
Frequency Array (LOFAR) (Rottgering 2003), the Square Kilometer Array (SKA) (van de 
Weygaert & van Albada 1996), the Mileura Widefield Array (MWA0) and the 21 cm Array 
(21CMA) (Pen et al. 2004), are ongoing or being planned. 

A major challenge to detect the primordial 21 cm signals is from the contamination of 
radio signals of the foreground. The brightness temperature fluctuations of the 21 cm signals 
would be seriously contaminated by high-redshift radio point sources, free-free emissions of 
the IGM, and the noises from the artificial radio interference in the VHF band (Di Matteo 
et al. 2002, 2004; Oh & Mack 2003). The power spectrum of the brightness temperature 
fluctuations of the redshifted 21 cm map may even be dominated by the fluctuations of 
the foreground (Di Matteo et al. 2004). Thus, we should search for useful features of the 
primordial 21 signals to distinguish it from the noisy foreground. For instance, the patchy 
structures of the 21 cm brightness temperature field are found to be highly non-Gaussian, 
which would be helpful to identify the 21 cm signals from the noisy foreground (He et al. 
2004, Mellema et al. 2006). Moreover, the correlation between the 21 cm signals and other 
emissions of the sources would also be useful for identifying high redshift sources. 

In this context, some studies focused on the 21 cm signals from individual UV ionizing 
sources in the reionization epoch. Their features may provide a direct identification of 
the ionized patches of the reionization (e.g., Tozzi et al. 2000). Very recently, the 21 cm 
emission and absorption around isolated ionizing sources has been calculated (Chuzhoy et al. 
2006; Cen 2006; Chen & Miralda-Escude 2006). However, their results, especially the time- 
dependence of the profile of the 21 cm emission and absorption, still do not quite agree with 
each other. This problem comes partially from the fact that they dropped the time derivative 
term of the radiative transfer equation, i.e., the speed of light is taken to be infinite. This 
approximation would be reasonable if the retardation effect is negligible. However, it is not 
the case for the 21 cm problem, for which we note that the size of the HII regions and the 
age of the sources used in these calculations have already been close to, or even violated, 
the retardation constraint r < ct. The retardation effect, or the effect of finite velocity of 
light, is not trivial. The time and space dependencies of the ionized and heated regions are 
substantially affected by the retardation effect (White et al. 2003; Shapiro et al. 2006a; Qiu 
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et al 2007). Therefore, it is necessary to re- investigate this problem by fully considering the 
kinetics of photons in the phase space. 

We will use the WENO algorithm for the radiative transfer calculation, which has been 
developed by Qiu et al. (2006, 2007). The WENO algorithm has been shown to have 
high order of accuracy and good convergence in capturing discontinuities and complicated 
structures (Shu 2003). The WENO algorithm is able to incorporate the temperature and 
species equations and to provide stable and robust solutions. It is also able to calculate the 
fractions of neutral hydrogen and helium within the ionized spheres, which generally are 
extremely small, but would be important to understand the possible correlation of the 21 
cm and Lja signals. 

In §2 we present the ionization and temperature profiles of the IGM around a UV 
ionizing source. §3 shows the 21 cm emission and absorption regions, and their dependency 
on intensities, frequency-spectrum and life time of the UV photon sources. Discussions and 
conclusions are presented in §4. The algorithm is summarized in the appendix. 



2. Ionization and Heating of a Point Source 
2.1. The model 

The whole process of reionization is very complex due to the inhomogeneous and clus- 
tering nature of the structure formation (e.g. Ricotti et al. 2002; Razoumov & Norman 
2002; Ciardi et al. 2003). Therefore, the configuration of the ionized regions generally would 
be non-spherical. However, in order to be comparable with previous studies (Chuzhoy et al. 
2006; Cen 2006; Chen & Miralda-Escude 2006), we will consider the ionization and heating 
of a uniform IGM around a point source. This simplified model is reasonable to reveal the 
features of the formation and evolution of the 21 cm emission and absorption regions. The 
radiative transfer equation of the specific intensity J(t, x, v) is (e.g. Bernstein 1988) 




where v is the photon frequency and H is the Hubble parameter at time t. Supposing the 
source is at the center r = \x\ = 0, we have S(t, \x\,v) = E(v)8(x), where E{y)dv is the 
energy of photons emitted from the central source per unit time within the frequency range 
from v to v + dv. The source term can equivalently be described by a boundary condition 

as 

lim47rr 2 J(t, r,u) = E{v). (2) 
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We will consider two kinds of frequency sepctrum: a power law spectrum E{y) = E (u /h / ) a 
with the total intensity (energy per unit time) E = E(y)dv = E u /(a — 1) and a 
blackbody spectrum. The effect of the expansion of the universe on the number density is 
negligible if the time scale concerned is a few Myr at the epoch 1 + z ~ 20. However, the 
effect of the redshift of the photon energy is essential to calculate the flux of Lya photons 
(Cen 2006). Therefore, we will keep the H term on the right hand side of eq.(l). 

The absorption coefficient in eq.(l) is given by k v (t,r) = a u (t,r)nu(t,r) and 

-(^X>M^, (3) 

where % runs over the components HI, Hel and Hell. <Ji{y) and rii(t,r) are, respectively, 
the cross section of ionization and the number density of component %. The number density 
rii is calculated by solving the species equations involving ionization and recombination 
processes. We use the on-the-spot approximation to take into account the effect of the 
diffusing photons and treat the coupling of H and He by the radiation of Hel recombination 
continuum and bound-bound transitions following the method in Osterbrock (1989). The 
corresponding coefficients for the ionization and recombination processes are taken from 
Fukugita & Kawasaki (1994) and Theuns et al. (1998). 

The IGM is heated by the UV photons with a rate in unit erg cm -3 s^ 1 as 

f°° v — v 
H = Y j n l duJ(t,r,u)a(u) ^(/ e ,i/) (4) 

i J »i v 

where v± is the threshold of frequency to ionize component %. The fraction of photon-electron 
energy converted to heat r)(f e , v) depends on the free electron fraction f e and the photon 
frequency v. We use the fitting formula of r](f e , v) given by Shull & van Steenberg (1985). 
The corresponding coefficients for radiative cooling processes are also taken from Theuns et 
al. (1998). 

The algorithm of the numerical solution of the radiative transfer equation coupled with 
the temperature and species equations of HI, Hel and Hell are described in Qiu et al. (2006, 
2007). A fifth-order finite difference WENO scheme and a third-order TVD Runge-Kutta 
time discretization are applied to the spatial and time derivative, respectively. A semi- 
implicit method is used for the temperature and species equations. The details and tests 
of the numerical algorithm are described in Qiu et al. (2006, 2007) and we give a brief 
description of the algorithm in the appendix. 
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2.2. Ionized sphere 

We first calculate the profile of an ionized sphere around a point source at 1 + z — 20. As 
an example, we consider a source of a power law spectrum with E = 7.25 x 10 44 erg s _1 and a 
spctral index a — 2, which are the typical values of QSOs at the EUV band (e.g. Zheng et al. 
1997). The initial ionization fraction of hydrogen and helium is set to be 10~ 4 , which is the 
remaining ionization fraction after recombination, and the initial temperature of the baryonic 
gas is set to be 10 K. A ACDM cosmology model with Q m = 0.27, Q b = 0.044, and h = 0.7 
is adopted. The profiles of the ionized sphere at time 0.575, 1.15, 1.725, 2.3 and 3 Myr are 
shown in Figure 1, in which / m = n m /n H , f HU = n HU /n H , / He i = n He i/n Ke , f ReU = n HeU /n Hc 
and /hcIii = n ncm /n Uc , where n H = n ffl + n H n and n Hc = n He i + n He ii + n He in. 

Figure 1 shows that although the main ionization zones have steep eges (I- front), there 
are still a small fraction of hydrogen and helium been ionized beyond the I-fronts. This 
ionization is produced by high energy photons which can effectively penetrate into the IGM 
further than photons near the ionization thresholds. The fraction /hii(^) has a long tail 
before it approaches the initial value 10 -4 ; the fraction of Hell shows a similar behavior as 
that of HII. The I-fronts of Hell and Helll generally are different from that of HII, and the 
decrease of the fraction of Helll with r is much faster than that of HII and Hell. 

We define the I-front rj(t) as the location where fm{ r i,t) = 10~ 4 . The speed of the 
I-front dri(t)/dt is plotted as the solid line in Figure 2. It shows that the evolution of the 
ionized sphere has two phases separated by a characteristic time t c . In the first phase t < t c , 
drj(t)/dt ~ c, and in the second phase t > t c , the speed decreases with t in a power law, 
drj(t)/dt oc t~ 3 / 4 . This feature of two phases was first pointed out by White et al. (2003) 
with an approximate solution by neglecting recombinations. This feature is actually not 
sensitive to the choice of the number 10~ 4 ; if we take fmij'ht) = 90%, the two phases 
feature still holds (Qiu et al. 2007). 

Within the I-front, / H i is on the order of 10~ 5 — 10~ 6 and corresponds to a column 
density of the order of 10 14 cm~ 2 on a physical length scale of the order of 100 kpc. This 
value of column density is just of the same order as the weak Lya absorptions (Bi & Davidsen 
1997). Therefore, the ionized sphere is not completely Gunn- Peterson transparent, but may 
produce Lya absorption features on the spectra of background sources. On the other hand, 
within the ionization front, fneii — 10~ 4 , and therefore, the Gunn-Peterson optical depth of 
Hell Lya is larger than that of HI. It leads to significant fluctuations of the ratio between 
the two optical depths (e.g. Shull et al. 2004; Liu et al. 2006). 

We also calculate the ionized sphere of sources with intensities E = 7.25 x 10 43 and 
7.25 x 10 45 erg s™ 1 and a power law index a = 2. Their I-front speeds are also plotted in 
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Figure 2. It shows that all the three curves of dri(t)/dt have the two phases. Especially, in 
the second phase the power law drj(t)/dt oc t~ 3 / 4 is approximately scaled with respect to 
(E) 1 / 2 , in other words, the characteristic time t c oc (E) 1 / 2 . This scaling can not be found 
with static solutions of the radiative transfer equation or by assuming the speed of light is 
infinite (e.g. Cen, ApJS, 2002; Whalen & Norman 2006). 

2.3. Kinetic temperature distribution 

Figure 3 shows the corresponding kinetic temperature distribution for the ionizing source 
with E = 7.25 x 10 44 erg s -1 and a power law index a = 2 at time 0.575, 1.15, 1.725, 2.3 
and 3 Myr. The kinetic temperature basically keeps constant around 3 x 10 4 K within the 
ionized sphere, which is determined by the balance between the photon-ionization heating 
and radiative cooling. An interesting feature is that the temperature at the center (r 0) is 
lower than the constant temperature. This is because the ionizing photon flux at the center 
is very strong, the number of /hi at r ~ is extremely low (see Figure 1), and therefore, the 
heating rate is low. 

Figure 3 shows that the gas is also substantially heated beyond the I-front. Actually, 
the energy spectrum of the radiation is significantly hardened around the I-front, because 
soft photons are effective to ionize the IGM, and hard photons are more effective to heat 
the IGM (Qiu et al 2007). After the time t c , soft photons are exhausted within the ionized 
sphere r < 77 (t); in contrast, hard photons are still abundant in the region r > rj(t), as hard 
photons can penetrate further beyond the ionized sphere. The region between the I-front 
and the front of light r = ct is a pre-heating layer. The temperature profile is different at 
different time and more elongated at a later time. This is because the speed of the I-front 
is much less than the speed of light after t c , while the spped of the front of heating is about 
the same as the speed of light. 

We introduce a radius r 50 (t) by the solution of T(r, t) = 50 K, which is about the CMB 
temperature at 1 + z = 20. We call the high temperature region to be rj(t) < r < r 50 (t) 
and the low temperature region r^o(t) < r < ct. The latter will play the role of the 21 
cm absorption. Since the speed of the radius r 50(f) is also less than c, the low temperature 
region increases quickly with time. Therefore, the kinetic temperature distribution shows a 
long tail in the range T < 50 K. We also calculate the kinetic temperature distribution of 
sources with intensities E = 7.25 x 10 43 and 7.25 x 10 45 erg s^ 1 and a power law index a — 2. 
The profiles of the temperature distribution is substantially ^-dependent. 
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3. 21 cm Emission and Absorption 

3.1. Spin temperature and 21 cm brightness temperature 

The 21 cm (z/ = 1420 MHz) line emission or absorption associated with the spin-flip 
transition of the electron in neutral hydrogen is significant if the number density n m of 
neutral hydrogen is high and the spin temperature T s , which describes the number ratio 
between the spin-up and the spin-down states, deviates from the temperature of cosmic 
microwave background (CMB) T cm ^. The preheating region discussed in §2.3, i.e., the region 
between the I-front and the front of light r = ct, would satisfy these two conditions. We 
may call this region as a 21cm shell. 

We do not consider the bulk motion of the baryonic fluid. Thus, when T s is much larger 
than the hydrogen hyperfine energy T* = 0.06 K, the optical depth of the 21 cm absorption 
r(z) is given by (Wild 1952; Field 1959) 



r(z) 



3hc 3 A 10 n K f m (z) 
32irvik B T s (z)H(z) 



(5) 



2.7 x 10' 3 f m (z) 



T cm b( z ) 



1/2 



Since t(z) <C 1, the observed brightness temperature excess (or deficit) for the mean gas 
density of the universe at the redshifted frequency v — vq/{1 + z) is 



5T b ~ 8 x 1(T 3 /hi 



T, 



cmb 



[l + z) 



1/2 



K. 



(6) 



The decoupling of the spin temperature T s from the CMB temperature T cmb could be 
achieved by the collision of neutral hydrogen and scattering of Lja photons, which is the 
Wouthuysen-Field mechanism (Wouthuysen 1952; Field 1959). The spin temperature T s is 
then given by a weighted average of the CMB temperature T cm b, the kinetic temperature T^, 
and the color temperature T c of background photons at Lya frequency (Field 1958, 1959) 

rp ^cmb VaT c ~\~ UcXk r-,\ 

s ~ i + y a + y c ' [) 

where y a is the Lya photon coupling coefficient and y c is the collision coupling coefficient. 

At redshift z = 20 — 30, the collision coupling is only efficient in overdense regions (e.g. 
Furlanetto & Loeb 2002; Shapiro et al. 2006b), and is negligible for cosmological mean 
densities. In the low density IGM, the Wouthuysen-Field mechanism would play the major 
role to decouple T s from T cmb . The resonant scattering of Lya photons by neutral hydrogen 
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leads to that the color temperature of the radiation spectrum near Lya frequency approaches 
the kinetic temperature of the baryonic gas. Therefore, if there are enough Lya photons, we 
have T s ~ T^. There are two important ways to contribute Lya photons: One is from the 
continuum spectrum photons redshifted to Lya frequency, which is governed by the H term 
in eq.(l). The other one is from the excitation of neutral hydrogen by high energy electrons, 
and we take the fitting formula of the fraction of excitation for high energy electrons from 
Shull & van Steenberg (1985). 



3.2. Evolution of the 21 cm shell 

Figure 4 shows the profiles of the spin temperature T s and the brightness temperature 
5T b for a point source with E = 7.25 x 10 44 erg s _1 and a power law index a = 2 at time 

0. 575, 1.15, 1.725, 2.3 and 3 Myr. The spin temperature T s for Tk > 10 4 K is set to be 
this will not affect our results since HI is highly ionized for T k > 10 4 K. In the 21 cm shell, 

1. e., the region between the I- front and the front of light r = ct, the brightness temperature 
5T b 7^ 0. In the early stage, when the speed of the I-front dri{t)jdt is close to the speed of 
light, there is no 21 cm shell. The 21 emission and absorption regions become broad when 
the difference between rj(t) and r = ct is significant. 

Figure 4 shows that the 21 cm layer always contains an emission region ST b > and an 
absorption region 5T b < 0. The maximum brightness temperature in the emission region is 

~ 35 niK, while the minimum temperature of the absorption region is ~ —100 150 mK. 

The absorption region grows faster than the emission region. Before 1 Myr, the emission 
region is comparable with the absorption region, but at 3 Myr, the emission region is about 
2 comoving Mpc while the absorption region exceeds 10 comoving Mpc. 

A prominent feature of the brightness temperature in Figure 4 is that the absorption 
region has a sharp cut-off at the front of light r = ct. This comes from the fact that in 
our model the Wouthuysen-Field coupling is efficient within the entire range of the photon 
propagation; consequently, in the region r^t) < r < ct, the Lya photons are abundant 
enough to lock T s with These photons come from the redshifted photons originally in 
the energy range between ionization energy and Lya line. In this case, the edge of the 21 
cm absorption region is an indicator of the lifetime of the ionzing source. 
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3.3. Dependence on the source intensity 

We also calculate the profiles of the spin temperature T s and the brightness temperature 
5Tb of sources with intensities E = 7.25 x 10 43 and 7.25 x 10 45 erg s -1 and a power law index 
a — 2. The results are presented in Figure 5. Comparing Figures 4 and 5, one can see 
that the maximum (emission) and minimum (absorption) of ST b basically are independent of 
the source intensity, while the profile of the 21 cm brightness temperature and its evolution 
depend strongly on the intensity. This is because the heating processes are sensitive to the 
evolution of photons in the phase-space, which also cannot be described by a static approach. 

Figure 6 shows the ratio between the sizes of emission and absorption regions, which 
decreases with time, because the growth of the absorption region (brightness temperature 
5T b < 0) is faster than that of the emission region. For a source with E = 7.25 x 10 45 erg 
s -1 , the size of the emission region is larger than that of the absorption region at an early 
stage t < 1 Myr, while for sources with E = 7.25 x 10 44 and 7.25 x 10 43 erg s -1 , the size of 
the emission region is always smaller than that of the absorption region. 



3.4. Dependence on the photon energy-spectrum 

Figure 7 shows the profiles of the spin temperature T s and the brightness temperature 
STj, for an ionizing source at 1 + z = 20 with the spectral index a = 3 and 4. The intensities 
for both cases are normalized to be E = 7.25 x 10 44 erg s _1 , i.e., the same as the source of 
§2.2. Comparing Figures 3 and 7, we see that the size of the 21 cm emission region of both 
a = 3 and 4 is narrower than that of a — 2. This is due to the lack of high energy photons 
for heating. On the other hand, the absorption regions are more notable than that in the 
case of a = 2, as the sources of a = 3 and 4 have stronger Lya photon fluxes than a — 2. 

The lights from first galaxies may have a thermal spectrum. Let us consider first galaxies 
consisting of stars with a blackbody spectrum and an effective temperature T ~ 10 5 K, which 
corresponds to population III stars with mass 200-500 M (Schaerer et al. 2002). Thermal 
radiation at temperature T ~ 10 5 K does not contain many photons with energy larger than 
100 ev. Therefore, the pre-heating beyond the I-front is very weak, and the 21 cm region 
is small. Figure 8 shows the spin temperature and the brightness temperature for a first 
galaxy with the same luminosity as the source for Figure 4. We see that the size of the 21 
cm emission region of Figure 8 is much narrower than that of Figure 4. However, the 21 cm 
absorption region is still significant, because the Lya photons are still efficient to support 
the Wouthuysen-Field mechanism. 
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3.5. Sources of short lifetime 

We now consider a source which is the same as that used for Figures 1 and 3, but the 
source stops to emit photons after 3 Myr. Figure 9 shows the kinetic temperature T^, the 
spin temperature T s and the brightness temperature <5TJ, at time 3, 3.45, 4.03 and 4.60 Myr. 
We see that the I-front and r 50 (t) still move ahead in the period of 3 to 3.45 Myr, but they 
slow down after 3.45 Myr. It is simply due to the fact that there are no more photons to 
support ionization and heating. Consequently, the 21 cm emission shell shrinks and almost 
disappears at 4.03 Myr. 

An interesting result is that the 21 absorption region is alive even when t = 4.6 Myr, 
i.e., 1.6 Myr after the photon source has died out. This is once again due to the abundance 
of Lya photons. When t > 4.03 Myr, no high energy photons can reach r > 6 Mpc, while 
the Wouthuysen-Field mechanism is working to the size r = ct. Figure 9 shows that the 
absorption region at t — 4.60 Myr is limited within two sharp cut-offs: T\ < r < r 2 . It is 
easy to check that T\ = c x (4.6 — 3) (Myr) and r 2 = c x (4.60) (Myr), which is just the range 
given by the retardation of the first and last light signals from the source with a luminous 
lifetime of 3 Myr. 



3.6. Thermal Broadening 

While including the thermal broadening, eq.(5) should be replaced by 

3hc 3 A 10 n m {z') 



rlz 



^ f Um[z ) F(z Adz' (8) 
^k B J T s (z>)H(z>) * {Z > Z) ^ 



where the factor F(z, z') is the normalized Doppler broadening line profile as 

with b = (2k B T '/mc 2 ) 1 ' 2 • For gas with temperature ~ 10 3 K, we have b ~ 1 x 10 -5 , and the 
comoving scale is thus D = 1 x 10~ 5 (1 + z)c/H(z). In the ACDM model at (1 + z) = 20, 
we have D ~ 0.02 h _1 Mpc. If the thickness of the 21 cm emission shell is less than D, the 
Doppler broadening will smooth out the 21 cm emission. As discussed in previous sections, 
the 21 cm emission regions are narrow for weak sources and high spectral index. Therefore, 
the 21 cm emission would not be observable for sources with E much less than 10 43 ergs s^ 1 
or E ~ io 43 ~ 45 ergs s" 1 , but a > 2. 
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4. Discussions and Conclusions 

We have studied the formation and evolution of the ionized and heated region around 
a UV ionizing source in the reionization epoch. In the first stage, the ionizing and heating 
fronts of the IGM are coincident and propagate with a speed close to the speed of light. The 
evolution enters the second stage when the frequency spectrum of UV photons is hardened 
due to the loss of photons by reionization. In this stage, the propagating speed of the I-front 
is less than the speed of light, but high energy photons and Lya photons can still reach the 
front of light r = ct. The spherical shell between the I-front rj(t) and r = ct is the region in 
which 21 cm signals are produced. The inner shell from rj{t) to r 50 (t) gives 21 cm emissions; 
the outer shell from r 50 (t) to r = ct produces 21 cm absorptions. 

Both the 21 cm emission and absorption regions are sensitively dependent on the in- 
tensity, frequency-spectrum and life-time of the ionizing source. If the intensity is much 
less than E ~ 10 43 ergs s _1 and the spectral index a is larger than 3, the emission region 
is very narrow, and the signal will probably be smoothed out by thermal broadening. For 
sources with intensity E ~ XO 44-46 ergs s^ 1 and a spectral index a — 2, the comoving size 
^5o(t) — r 'i(f) is about 1-5 Mpc. The brightness temperature excess upon the CMB is about 
30 mK. Therefore, it yields a radio signal near the band of 70 MHz and is resolvable with 
observations of angular resolution better than 0.5 arcmin and spectral resolution better than 
40 kHz. 

The 21 cm absorption region has very different behaviors from the emission regions. 
Absorption region mainly depends on the coupling of the spin temperature and the kinetic 
temperature, which is efficient if there are enough Lya photons. The ionizing source with a 
power law spectrum or thermal spectrum at T ~ 10 5 can always provide enough Lya photons 
until the largest distance r = ct. That is, the size of the absorption region directly measures 
the age of the source. For a source with lifetime t s , the absorption region can survive even 
when t is larger than t s . In this region, the brightness temperature deficit with respect to 
the CMB is of the order of a few ten mK. The comoving size of the absorption region can be 
as large as 10 - 20 Mpc. Moreover, the small fraction of HI and Hell remained in the ionized 
sphere can lead to absorption features on the spectra of background sources. Therefore, the 
correlated 21 cm emission, absorption and Lya signals would be helpful to identify ionized 
patches in the early universe. 

Although all the calculations of this paper are on spherical regions, many features re- 
vealed with the spherical profiles would also hold for non-spherical profiles. This is because 
many features of the formation and evolution of 21 cm signals depend only on the growth of 
the ionized and heated regions, but not on the details of their configuration. For instance, 
it can be shown that the growth of a non-spherical ionized region around clustered sources 
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also consists of two phases as that shown in Figure 2. In the first phase, the speed of the 
I-front is also close to the speed of light, even though in these cases we cannot use one ionized 
radius 77 to describe the ionized region. Moreover, similar to the spherical ionized regions, 
the characteristic time t c is longer for clustered sources with higher intensity E. Therefore, 
the stronger the sources, the later the formation of the 21 cm regions. This property is 
regardless of the geometry of the clustered sources (Qiu et al. 2007, in preparation). 

In our calculations, the dynamical behavior of cosmic baryonic gas is ignored. Such an 
approximation would be reasonable when the speed of the sound wave are very different from 
the speed of the formation and evolution of the 21 cm emission regions. Most features of the 
21 cm emissions and absorptions should still hold when the hydrodynamics of the baryonic 
gas is taken into account. 
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A. Numerical Algorithm 

A.l. Approximation to the spatial derivative 

To approximate the spatial derivative in equation (1), the fifth-order finite difference 
WENO scheme is used. Specifically, to calculate dJ/dr, the variable v is fixed and the 
approximation is performed along the spatial coordinate r 

— J(t n , ri ,u) « — (h i+1/2 - hi- 1/2 ) , (Al) 

where the numerical flux /ij+1/2 on a uniform mesh r« at time step n is obtained with the 
procedure given below. We can use the upwind fluxes without flux splitting in the fifth-order 
WENO approximation because the wind direction is fixed (positive). 

First, we denote 

ht = J(t n ,r h is), i = -2,-l,...,N r + 2. 
The numerical flux from the regular WENO procedure is obtained by 

hi+l/2 = + ^2^+1/2 + w 3^S+i/2> 



where h^^ 2 are the three third order fluxes on three different stencils given by 

n i+i/2 ~ 3^-2 ~ Q n i-i + 



n i+i/2 ~ ~Q tli ~ 1 + 6 3 



and the nonlinear weights uo m are given by 
with the linear weights 7; given by 

1 _ 3 3 

7i - TJy 72 - g, 7s - Try 

and the smoothness indicators given by 

A = ^ (^- 2 " 2/l *" 1 + + \ (hi-2 - 4/ii-i + 3^) 2 

13 2 1 2 

^ 3 = 12 ^ _ 2ki+1 + ^ +2 ^ + 4 ^ ~~ 4/li+1 + hi+2 ' ' 
e is a parameter to avoid the denominator to become and is taken as e = 10 
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A. 2. Numerical calculation in v 

The computational domain in v is on a non-uniform mesh 

Vj = wif/i & = jAf , A£ = \og 2 u max /N v , j = 0, N u . 

The integration of heating rate eq. (4) and photon-ionization rate is approximated by a 
fourth order quadrature formula 

POD °° 

/ f(x)dx = Ax^wjfUAx) + 0(Ax 4 ), (A2) 



"° J=J0 



where z/o is the threshold frequency of ionization, and the weights Wj are given by 
3 7 23 

w jo = g, w io+l = g, w io+2 = Wjo+j = 1, for J > 2 - 
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A.3. Time derivative 

The time derivative dJ/dt is calculated using the third-order TVD Runge-Kutta dis- 
cretization: 

= u n + AtL(u n ,t n ) (A3) 
w (2) = ^ u n + ^(u {1) + AtL(u w )) (A4) 

u n+1 = -u n + -(u^ + AtL(u^)). (A5) 
3 3 

We adapt the multi-timescale strategy and use a semi-implicit scheme in the evolution of 
temperature and species. It greatly releases the time step restriction, therefore saves compu- 
tational cost. We refer to Qiu et al. (2007) for more details of the numerical implementation. 
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Fig. 1.— /hi, /hii, /Hei, /Heii, and /hciii, the fractions of HI, HII, Hel, Hell, Helll around a 
point source at redshift 1 + z = 20 with an intensity E = 7.25 x 10 44 erg s -1 and a power low 
index a = 2. From left to right, the time is 0.575, 1.15, 1.725, 2.3 and 3 Myr, respectively. 
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Fig. 2. — The speed of the I-front dri(t)/dt for three different intensities E. 
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Fig. 3. — The profiles of the kinetic temperature T(t,r) for the same model as in Figure 1. 
From left to right, the time is 0.575, 1.15, 1.725, 2.3 and 3 Myr, respectively. 
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Fig. 4. — The profiles of the spin temperature T s and the brightness temperature ST b for the 
same model as in Figure 1. From left to right, the time is 0.575, 1.15, 1.725, 2.3 and 3 Myr, 
respectively. 
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Fig. 5. — The profiles of the spin temperature T s and the brightness temperature 5Tb for 
point sources with a spectral index a = 2, but with intensities E = 7.25 x 10 43 (left) and 
7.25 x 10 45 erg s" 1 (right). In each figure, from left to right, the time is 0.575, 1.15, 1.725, 
2.3 and 3 Myr, respectively. 
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Fig. 6. — The time dependence of the ratio between the size of emission and absorption 
regions for three different intensities with spectral index a — 2. 



-23- 





5 10 15 

r (comoving Mpc) 



20 



-50 



o -100 

-150 
-200 



urn 




5 10 15 

r (comoving Mpc) 



20 



Fig. 7. — The profiles of the spin temperature T s and the brightness temperature 5Tb for a 
point source with spectral indexes a = 3 (left) and 4 (right). The intensities are normalized 
to have the same number of ionizing photons as the model in §2.2. In each panel, from left 
to right, the time is 0.575, 1.15, 1.725, 2.3 and 3 Myr, respectively. 
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Fig. 8. — The profiles of the spin temperature T s and the brightness temperature for a 
source with a thermal spectrum of T = 10 5 K with the same luminosity as in Figure 4. In 
each panel, from left to right, the time is 0.575, 1.15, 1.725, 2.3 and 3 Myr, respectively 
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Fig. 9. — The profiles of the kinetic temperature Tk, the spin temperature T s , and the 
brightness temperature STj, for the same model as in Figure 4, but it stops to emit photons 
after t — 3 Myr. In each panel, from left to right, the time is 3, 3.45. 4.03 and 4.60 Myr, 
respectively. r\ and r 2 are the boundaries of the profile ST b at time 4.60 Myr. 



